Search for: All records

High-resolution vibrational spectra of C–H, O–H, and N–H stretches depend on both molecular conformation and environment as well as provide a window into the frequencies of many other vibrational degrees of freedom as a result of mode mixing. We review current theoretical strategies that are being deployed to both aid and guide the analysis of the data that are encoded in these spectra. The goal is to enhance the power of vibrational spectroscopy as a tool for probing conformational preferences, hydrogen bonding effects away from equilibrium, and energy flow pathways. Recent years have seen an explosion of new methods and strategies for solving the nuclear Schrödinger equation. Rather than attempt a comprehensive review, this work highlights specific molecular systems that we have chosen as representing bonding motifs that are important to chemistry and biology. We focus on the choices theoretical chemists make regarding the level of electronic structure theory, the representation of the potential energy surface, the selection of coordinates, preferences in basis sets, and methods of solution. 

The pathways and timescales of vibrational energy flow in nitromethane are investigated in both gas and condensed phases using classical molecular mechanics, with a particular focus on relaxation in liquid water. We monitor the flow of excess energy deposited in vibrational modes of nitromethane into the surrounding solvent. A marked energy flux anisotropy is found when nitromethane is immersed in liquid water, with a preferential flow to those water molecules in contact to the nitro group. The factors that permit such anisotropic energy relaxation are discussed, along with the potential implications on the molecule’s non-equilibrium dynamics. In addition, the energy flux analysis allows us to identify the solvent motions responsible for the uptake of solute energy, confirming the crucial role of water librations. Finally, we also show that no anisotropic vibrational energy relaxation occurs when nitromethane is surrounded by argon gas. 

Phenol–benzimidazole and phenol–pyridine proton-coupled electron transfer (PCET) dyad systems are computationally investigated to resolve the origins of the asymmetrically broadened H-bonded OH stretch transitions that have been previously reported using cryogenic ion vibrational spectroscopy in the ground electronic state. Two-dimensional (2D) potentials describing the strongly shared H atom are predicted to be very shallow along the H atom transfer coordinate, enabling dislocation of the H atom between the donor and acceptor groups upon excitation of the OH vibrational modes. These soft H atom potentials result in strong coupling between the OH modes, which exhibit significant bend-stretch mixing, and a large number of normal mode coordinates. Vibrational spectra are calculated using a Hamiltonian that linearly and quadratically couples the H atom potentials to over two dozen of the most strongly coupled normal modes treated at the harmonic level. The calculated vibrational spectra qualitatively reproduce the asymmetric shape and breadth of the experimentally observed bands in the 2300–3000 cm–1 range. Interestingly, these transitions fall well above the predicted OH stretch fundamentals, which are computed to be surprisingly red-shifted (<2000 cm–1). Time-dependent calculations predict rapid (<100 fs) relaxation of the excited OH modes and instant response from the lower-frequency normal modes, corroborating the strong coupling predicted by the model Hamiltonian. The results highlight a unique broadening mechanism and complicated anharmonic effects present within these biologically relevant PCET model systems. 

We report and interpret recently recorded high-resolution infrared spectra for the fundamentals of the CH 2 scissors and CH stretches of gas phase cyclopentane at −26.1 and −50 °C, respectively. We extend previous theoretical studies of this molecule, which is known to undergo barrierless pseudorotation due to ring puckering, by constructing local mode Hamiltonians of the stretching and scissor vibrations for which the frequencies, couplings, and linear dipoles are calculated as functions of the pseudorotation angle using B3LYP/6-311++(d,p) and MP2/cc-pVTZ levels of theory. Symmetrization ( D 5 h ) of the vibrational basis sets leads to simple vibration/pseudorotation Hamiltonians whose solutions lead to good agreement with the experiment at medium resolution, but which miss interesting line fractionation when compared to the high-resolution spectra. In contrast to the scissor motion, pseudorotation leads to significant state mixing of the CH stretches, which themselves are Fermi coupled to the scissor overtones. 

Correction for ‘Single-conformation spectroscopy of cold, protonated D PG-containing peptides: switching β-turn types and formation of a sequential type II/II′ double β-turn’ by John T. Lawler et al. , Phys. Chem. Chem. Phys. , 2022, 24 , 2095–2109, https://doi.org/10.1039/D1CP04852J. 

d -Proline ( D Pro, D P) is widely utilized to form β-hairpin loops in engineered peptides that would otherwise be unstructured, most often as part of a D PG sub-unit that forms a β-turn. To observe whether D PG facilitated this effect in short protonated peptides, conformation specific IR–UV double resonance photofragment spectra of the cold (∼10 K) protonated D P and L P diastereomers of the pentapeptide YAPGA was carried out in the hydride stretch (2800–3700 cm −1 ) and amide I/II (1400–1800 cm −1 ) regions. A model localized Hamiltonian was developed to better describe the 1600–1800 cm −1 region commonly associated with the amide I vibrations. The CO stretch fundamentals experience extensive mixing with the N–H bending fundamentals of the NH 3 + group in these protonated peptides. The model Hamiltonian accounts for experiment in quantitative detail. In the D P diastereomer, all the population is funneled into a single conformer which presented as a type II β-turn with A and D P in the i + 1 and i + 2 positions, respectively. This structure was not the anticipated type II′ β-turn across D PG that we had hypothesized based on solution-phase propensities. Analysis of the conformational energy landscape shows that both steric and charge-induced effects play a role in the preferred formation of the type II β-turn. In contrast, the L P isomer forms three conformations with very different structures, none of which were type II/II′ β-turns, confirming that L PG is not a β-turn former. Finally, single-conformation spectroscopy was also carried out on the extended peptide [YAA D PGAAA + H] + to determine whether moving the protonated N-terminus further from D PG would lead to β-hairpin formation. Despite funneling its entire population into a single peptide backbone structure, the assigned structure is not a β-hairpin, but a concatenated type II/type II′ double β-turn that displaces the peptide backbone laterally by about 7.5 Å, but leaves the backbone oriented in its original direction.